Annihilation catastrophe: Prom formation to universal explosion 
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I present a systematic tlieory of formation of the universal annihilation catastrophe which develops 
in an open system where species A and B diffuse from the bulk of restricted medium and die on its 
surface (desorb) by the reaction A + B —i- 0. This phenomenon arises in the diffusion-controlled limit 
as a result of self-organizing explosive growth (drop) of the surface concentrations of, respectively, 
slow and fast particles {concentration explosion) and manifests itself in the form of an abrupt singular 
jump of the desorption flux relaxation rate. As striking results I find the dependences of time and 
amplitude of the catastrophe on the initial particle number, and answer the basic questions of when 
and how universality is achieved. 

PACS numbers: 82.20.-w, 05.70.Jk, 05.70.Ln 



I. INTRODUCTION 

For the last decades it has been shown that in spite 
of the fundamental simplicity the reaction-diffusion sys- 
tem A + S — > exhibits rich cooperative behavior 
0. One of the most impressive examples in the class 
of systems with bulk reaction is the phenomenon of 
Ovchinnikov-Zeldovich segregation (spontaneous growth 
of single species domains which leads to anomalous reac- 
tion deceleration). Here I focus on another wide class of 
systems where reaction proceeds on the catalytic surface 
of medium whereas diffusion proceeds in its bulk. In the 
work 2| it was first demonstrated that in this class of 
systems the interplay between reaction and diffusion ac- 
quires qualitatively new features and leads to a new type 
of self-organization. It has been found that once particles 
A and B diffuse at different mobilities from the bulk of 
restricted medium onto the surface and die on it (desorb) 
by the reaction A + B ^ 0, there exists some threshold 
difference in the initial numbers of A and B particles, Ac, 
above which the loop of positive feedback is "switched 
on" and the process of their death, instead of the usual 
deceleration, starts to accelerate autocatalytically. Re- 
cently, it has been discovered Q that the deceleration- 
acceleration transition is a prelude to much more non- 
trivial dynamical effects : in the diffusion-controlled limit 
A — > cx) a new critical phenomenon develops - annihila- 
tion catastrophe, which arises as a result of self-organizing 
explosive growth (drop) of the surface concentrations of, 
respectively, slow and fast particles [concentration ex- 
plosion) and manifests itself in the form of an abrupt 
singular jump of the desorption flux relaxation rate. 

The key features of the annihilation catastrophe have 
been obtained on the assumption that the initial num- 
ber oi A — B pairs is large, so that after a transient 
stage the annihilation dynamics crosses over to the in- 
dependent of initial pair number universal regime 0, 
0. Until now, however, the principal questions remain 
open: When and how is the universal regime achieved? 
Moreover, the central question concerned with the time 
moment of the catastrophe as a function of the initial 
pair number remains open too. In this paper I present a 



systematic theory which gives the exhaustive answers to 
these questions. 



II. MODEL 

I consider a model, in which species A and B are sup- 
posed to be initially uniformly distributed in the bulk 
of an infinitely extended slab of thickness 2i. Both 
species diffuse to the surface X = ±£ (X e [--(., P\) 
and irreversibly desorb as a result of surface reaction 
Aads + Bads AB with the rate proportional to the 
product of surface concentrations / = kcasCbs 0(Fig. 
1). Because of planar spatial homogeneity the system 
is effectively one dimensional. The boundary conditions 
are determined from the equality of diffusion and des- 
orption / flux densities at the surface l'^\s — I, i.e., it is 
assumed that the surface layer capacity can be neglected. 
According to 0, after introducing index "H" (heavy) 
for slower diffusing species and index ^^L" (light) for a 
faster one, the problem of species evolution in the dimen- 
sionless units reads (by symmetry I consider the interval 
[0,£] only) 



dh/dr ■ 



dl/dr = (l/p)V^l, (la) 
(l/p)Vl ~hj,, (lb) 



with V(h,l) \x=o~ and initial conditions h(x,0) = Hq 
and l(x,0) ~ Iq. Here h(x,T) — ch/ci, and 1(x,t) = 
cl/c^ are the reduced concentrations, V = d/dx,x — 
X/£ E [0, 1] is the dimensionless coordinate, r — Dnt/i^ 
is the dimensionless time, p = Dh/Dl < 1 is the ratio of 
diffusivities, and c* = Dh/k£ is the characteristic con- 
centration scale. According to (lb) particles disappear 
in pairs only, i.e., 

J = hJs = -< h> ^ -< I >, 

where J — I/I^, is the reduced desorption flux density 
and /* = Kcl is its characteristic scale, therefore 



< h> - < l >= A 



const. 
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i.e., the excess amount stays "inert" in the bulk (here 

< h>= hdx = Mh/M* and < I >= Idx = Ml/M* 
are the total reduced numbers of particles in the bulk per 
unit of surface and AC = = Dh/k). This "inert" part 
of the majority species A = SAf/Af», acts as a control pa- 
rameter, whereas its "active" part A'' = Afpair/-^* (equal 
to the total mimber oi H — L pairs) acts as the only 
variable decaying from A'o to as r ^ oo. According 
to [2,3], the key features of nontrivial dynamics, devel- 
oping in the system (1) in the diffusion-controlled limit 
-^0 = cl{0)k£/Dh 00, may be formulated as follows: 
When parameter A achieves the threshold value 

Ac = \/uiQ/pta.n{^/w^) (2) 

(wo = 77^/4 is the main eigenfrequency of the diffusion 
field relaxation) , the system undergoes a transition to the 
state where after a transient stage the surface concentra- 
tion of i?-particles, hg, and as a result the rate of death 
of pairs start to grow with time autocatalytically (growth 
of hs accelerates the drop of Is, the drop of Is accelerates 
the growth of hs, and so on). With growing A the au- 
tocatalytic stage becomes more and more pronounced so 
that far beyond the threshold self-acceleration acquires 
explosive character: at A ^ oo the rates of growth 
V,Hs = +^^57^ and relaxation Qls = -^^7^ are syn- 
chronized singularly growing by the law 

ns = l/\T\, |r| = |r-r.|^0, 

where the point of finite-time singularity r^, is achieved 
at the moment when the reduced number of pairs n(r) = 
N{t)/ A drops to some critical value n^,. The most spec- 
tacular consequence of concentration explosion is singu- 
lar behavior of the flux relaxation rate tJ^ = — ^^r^ 
which is siistained constant up to the critical point t^, 
upon reaching which tJ^ blows up abruptly to 00: at 
K. = p^/^A/Ac — > 00 the width of the jump is contracted 
and its amplitude grows by the laws 

|T|cat cx A-2/5 ^ 0, maxrj^ cx AV* ^ 00. 

In the work [3] it has been shown that in the limit 

no = iVo/A 00 this catastrophical jump of tJ^, called 
annihilation catastrophe, proceeds in the universal (no - 
independent) regime and the scaling theory of universal 
explosion has been given. However, the approach de- 
veloped in [3] did not allow one to say anything neither 
about the dynamics of explosion formation, nor about 
how the universal regime is achieved, nor about how the 
point of catastrophe depends on the initial conditions. 
The goal of this paper is to give a closed theory of an- 
nihilation catastrophe formation and, based on it, (a) to 
reveal the catastrophe universalization regularities and 
(b) to find the dependence Ti,(no) for arbitrary ratio of 
diffusivities. I show that this strongly nonlinear prob- 
lem allows for strict and elegant analytical solution, I re- 
veal its surprisingly rich " structure" , and demonstrate 
remarkable agreement with numerical calculation results. 



III. THEORY OF UNIVERSAL ANNIHILATION 
CATASTROPHE FORMATION 

I start with the exact formal solution of the problem 
(1) in the Laplace space f{s) = Cf{T) 

h(x,s) = hQ/s + (hs — ft-o/s) cosh(x-y/s)/ cosh(\/fl). 
l{x, s) = Iq/s + {Is — lo/s) cosh.{x^/sp)/ cosh(y^). (3) 

According to (3) the BC's (lb) acquire the form 

J = {ho /s — hs)^/s tanh \/s = 
= {lo/s - fs ) v^sTp tanh = C.{hsls) (4) 

and in an implicit form completely define the behavior of 
surface concentrations hs and Is which in turn via Eqs.(3) 
define the evolution of spatial distribution. The strategy 
for solution of the nonlinear chain (4) is based on that in 
the ff-diffusion-controUed regime the ratio should 
rapidly drop with the time, therefore, according to (4) 
we can first (i) calculate J(r) and Is{t) neglecting the 
hg contribution, then (ii) derive ^.^(t) from the condi- 
tion hs = J /Is, and finally (iii) calculate next- to- leading 
terms thereby defining the self-consistent picture of sur- 
face concentrations evolution. 



A. Transient dynamics r <C 1. 

At sufficiently small times the flux density is slightly 
changed so assuming J « Jq = holo from (4) one obtains 

hs = ho{l-Vh + ---), ls = lo{l-vi + ---) (5) 

where Vi = -^^r/n, th = l/ll and t; = I /phi- Ac- 
cording to (5) relative drop rate for Is and hs is governed 
by the value of the parameter 

R = vi/vh = r^/p 

where r = ho/lo = (1 + no)/no- So, the necessary condi- 
tions for the iJ-diffusion-controUed annihilation regime 
are fulfilment of the requirements ^o = -^0 S> 1 and 
R < Rc ~ 1- Taking the both requirements fulfilled 
from (5) and (4) one concludes that at r/( <C r <C 1 the 
flux should drop by the law 

J = ho{l-m)/^^ Jo{Th/Tf'W^ (6) 

with m \T/Th^oo~^ 0. In the region Th <^ t <^ p where the 
both species diffuse in the semi-infinite medium regime, 

from (4) it follows Ig — hs^ = lo^ where e = Rc — R. 
Assuming e > to be not too small from (4) and (6) we 
find 

hs = —^{l -(!> + ■■■), Is = loe{l + (!> + ■■■) (7) 

where (f) = {R/ ^€^)^Th/T. Using then (5) and (7) we 
self-consistently find m ~ (i?/e^)^(r/i/T) <C </> to conclude 
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that asymptotics (7) is realized at the condition e ^ 

T^^^ = \/-/Nq. In the opposite hmit < e <C r^^'* by 
the same procedure we come to the critical asymptotics 

-<? + •••) - VP^s(l + .g + ---) = 

where g ^ eijlTy^^l^ and — = ciry^jT^^I'' with 

c = 7ri/'*r(|)/2r(i). At sufficiently smah p and not too 
large r (so that i? <C i?c) in the region r;i,p ^ t <C 
l/r^ < 1 the L particles distribution becomes uniform. 
In this regime from (4) and (6) we find 

/i, = + a + -?o(l-<T + ---), (8) 

where u ^ 2r-\/ t/tt + 0{ry/Th, Ry^p/r). 

B. Self-accelerating dynamics t > 1. 

According to Eqs.(7), (8) at large Nq ^ oo, e ^ 
I/1//V0 and not too large r (i.e not too small no) by 
the moment r ~ 1 when the diffusion length of H parti- 
cles becomes comparable with the system's size, the ratio 
hg/ho oc r/cNo 0. Neglecting the hg contribution, it 
can be shown (see below) that at r > 1 and large no the 
hs value has to exponentially rapidly tend to a constant 
C. In view of this, according to (4) we write 

J ^ {ho- C) tanh(\/s)/Vs 
whence it follows 

J = yle"'^°^(l + e-*'^°^ + ---), (9) 

where A = 2{ho — C) w 2hQ. With the same accuracy 
from (4) we have 

is = Iq/s - [{ho - C)/s]^tanh(v^) coth(^/ip) 

whence it follows 

= (^/A,)e--"^(l-A), (10) 

where Ac is defined by Eq.(2) and the leading contribu- 
tion to A is governed by the sum of exponentially decay- 
ing, A_, and exponentially growing, A+, terms 

A = A_ + A+, 

A_ = gse-^'^«^ + Qpe-'^'^''\ A+ Ae"°^ (11) 
with exponent x— i^/p ^ 1) ^-nd amplitudes 
gs = ~{Ac/3)\/p/uJo cot (3^aJ^) , 



A = Ac(A-C)/A 

From the condition hs ~ J /Is and Eqs.(9),(10) locking 
the chain we find 

hs^A,{l + e-'^'>- + ■■■)/{! -A). (12) 

According to (11), (12) in the limit of large > 00 

(|A| 0) the hs value at any A exponentially rapidly 
achieves universal asymptotics h^ = Ac whence it fol- 
lows C = Ac- Essentially that at p < pc = 4/9 in A_ 
the first term dominates, therefore the relaxation rate 
is independent of p, whereas at p > pc in A_ the sec- 
ond term dominates, therefore the relaxation rate decays 
with growing p. In view of (7), (8), and (12) using the 
time convolution it is easy to check that the contribution 
of transient stage is reduced only to a relative shift of the 
amplitude 

6Atr/A Iho^oo-- {r/eho) / e'^"' dO ^ 
Jo 

therefore with an accuracy of vanishingly small terms we 
finally have 

A = Ac(A- Ac)/2/io. 

Two important consequences immediately follow from 
Eqs. (11), (12): (i) at A = Ac the amplitude A re- 
verses its sign which is the first rigorous proof of the 
threshold arising of self-acceleration for arbitrary ratio 
of diffusivities < p < 1 (note that at the long-time 
tail of self- acceleration r — s- 00, /i^ — s- A a strict deriva- 
tion of the threshold Ac has been given in the work 5]); 
(ii) for the time of self-acceleration start, Ts, and depar- 
ture of the starting point /i™'" from the critical asymp- 
totics hs, 6s = (/j™'" — Ac)/Ac, we find, respectively, 
Ts = [l/wo(a -I- 1)] ln(Q!g/A) and 

= (a+l)e^(A/a)^, (13) 

where in the p < pc range a ~ 8, g ~ 1 + gs whereas in 
the Pc < p < 1 range a = g — gp. 

C. Annihilation catastrophe. 

According to (4) and (12) the exponential growth 
Shs/hs = A+ leads to the exponentially growing con- 
tribution to the flux 5J/J = — /3A^ where (3 = 
-^oJo tanh Y^/(A — Ac). Far beyond the threshold /3 oc 
A^^ suggesting that in the limit A — > cxd at the ini- 
tial stage of self- acceleration the contribution to the fiux 
is vanishingly small. The remarkable fact to be proved 
below is that at A — > 00 the contribution to the fiux 
remains vanishingly small up to the point of finite-time 
singularity A(r*) — > 1 where hs/hs 00. This means 
that Eqs.(9)-(12) give complete description of the explo- 
sion dynamics. Moreover, as at A — > cxd the parameter 



Qp = (AcVp/7r)tan(7r/^), 



A = Ac/2(l-|-no) 
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becomes the unique function of no = Nq/A, Eqs.(9)-(12) 
allow one to achieve two main goals: a) to find the time of 
the catastrophe T^,{no) and b) to obtain the description of 
explosion evolution with growing hq. Taking A(t^,) = 1 
and A < 1 from Eqs.(ll) wc find 

r. = r:(l + <50, (14) 

where 

= (4/7r2) ln[2(l + no)/A,] 

and 

^r(no)--(feA« + f?pAX)/|lnA|. 

Introducing then a relative time T = t — t-i, from Eq.(12) 
we find that at any p < 1 in the vicinity |T| ^ uJq^ an 
explosive growth of hg sets in by the law 

K = il + Q)/^i\T\, |T|^0, (15) 

where fi = ^ ^ 1 — p and 

Q{no) = {l + 9gs)X^ + {l+x)e,A''- 

According to Eq.(9) in the course of explosion the flux is 
actually "frozen" 

J = J,[l + cjo(l + w)|T| + ---] 
reaching at the point of singularity the value 

J, = AA,(1 + G), |T|^0, (16) 

where u'(no) = 8A^ and 

G{no) = (l + ^?8)A« + ^^pA>^. 

From Eqs.(15),(16) and the condition hgls = J* there 
immediately follows the synchronization of the growth, 
^Hs = +^"^§7^, and relaxation, JIls = -^17^: rates 
which singularly grow by the law 

Qs = ^Hs = ^Ls = l/\r\. 

Clearly the dominant contribution to the explosion- 
triggered "antiflux" Jex occurs in the vicinity |T| <C 1 
where the diff'usional response to the explosion forms in a 
narrow layer cx ^/\T\ 0. Thus, considering the medium 
to be a semi-infinite one and allowing for (15) we can 
write 

_ _ dhs de (1 + Q) 

whence there follows smallness of \ Jex\/ J* down to |T| cx 
A-2/3 ^ 0. Calculating then a singular contribution into 
the flux relaxation rate 

^7'=-^=^o(l + u;) + [r7^] (17) 



we find 

whence there follows the catastrophic jump of tJ^ from 
tJ^ k ujq to tJ^ -^00 with the width 

According to the culmination consequence of the 
explosion in the limit K. ~ p^^^A/Ac — > 00 is the ex- 
act scaling description of passing through the point of 
singularity based on two key conditions: a) requirement 

^/phs = Is 

providing for equality of diffusional responses = J^^ 
(here and in what follows () = dQ/dr) and b) self- 
consistent condition of passing through the singularity 
point at a "frozen" flux 

It has to be mentioned, however, that the scaling be- 
havior of fig was postulated in [|| on a basis of indirect 
arguments and only later the postulated scaling function 
was analytically substantiated. Moreover, the space of 
being concise, an important chain of considerations 
remained beyond the analysis given there. Below, I shall 
give for the first time a systematic scaling theory of the 
annihilation catastrophe in the universal limit uq ~^ 00, 
and then, on its basis, a complete picture of catastrophe 
universalization with a growth of hq will be constructed. 



D. Scaling laws of concentration explosion and 
annihilation catastrophe. 

For simplicity, I shall begin with inferring the scal- 
ing laws of concentration explosion and annihilation 
catastrophe in the universal limit no — s- 00. Taking 
X,Q,G ^ 0, according to (9) and (10) at the explosion 
stage tiJolTj — + we have 

- j.(i+c.om + ...) 

and 

/w = M./*|T|(i + ^o|r|/2-f •••), 

whence it follows 

/.r = ^^"V^i"^ = i/^m + • • • , 

where the index "(0)" marks the solutions which neglect 

the contribution of hg = h'^^ (/li"'' =0). As it has been 
mentioned above, an explosive growth of h^^ ought to 
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trigger an explosive growth of "antiflux" J^x in calcu- 
lating which the medium can be regarded to be a semi- 
infinite one 



f 



de 



(19) 



(apparently, to an accuracy of vanishingly small terms 
the lower limit of the integral can be directed to — oo). 
Substituting here hi", we find 

Je. = -aj/^\T\'/^, 
where aj = So, for the total flux we have 



J = j(0) + j^^ = J Jl+a;o|T|- 



a.j 



MJ*|7'|3/2 



As the diffusion fluxes of fast and slow particles must 
be equal j£\s = J^ls = J, it is clear that against the 
background of dropping li^^ there must arise an explo- 
sive growth of which must initiate exactly the same 
"antiflux" = = Jex- Assuming that at a devel- 
oped explosion stage ClsP 3> 1 the dominant contribution 
to Jex occurs at times \T\/p -C 1, when for L particles 
the medium can be regarded to be semi-infinite, by full 
analogy with Eq. (19) we can write 



J p.: 



f 

J — OO 



dV:^ 



dO 



de ^p-k{t - e) 



(20) 



Thus, by requiring the equality of (19) and (20) we have 

If = vphT = vp/(^\n 

Making use of this result we obtain 



and finally from the condition hs = J /Is we find 



(21) 



(22) 



where /,/ = ajii^\r\/p. From Eqs. (21) and (22) it fol- 
lows that in the vicinity of some characteristic time 7} ~ 
]?- 1 I [i\f7l the explosion rate growth begins to drasti- 
cally decelerate. Due to the requirement << 1 
necessary for the realization of the synchronous explo- 
sion regime (20), the explosion deceleration begins long 
before a noticeable flux departure from the critical one 
\Jex^S^I J-k\ ^ 1- Introducing the parameter 

it can easily be seen that at any finite < p < 1 with a 
growth of A in the limit of large /C p^/^A/Ac — > oo 



|Je.(r/)|/J.~ JT//P~1//C 



-1/4 



0. 



and, therefore, down to 7}' ^ the flux remains "frozen" . 
One of the most important consequences of the drastic 
deceleration of the explosion rate growth is the drastic 
deceleration of the flux relaxation rate growth. As it will 
be shown in what follows, in the limit of large /C — » oo 
as a result of such deceleration the flux remains "frozen" 
and, therefore, the explosion develops synchronously 



t-Ls 



(23) 



hoih before and after passing through the point of sin- 
gularity where fig reaches a maximum. I shall now show 
that the condition (23) along with the following from (19) 
and (20) key condition 



It 



VpK 



(24) 



lead to a remarkably complete scaling description of pass- 
ing through the point of singularity. Taking 1^ = li^^ +1^^ 
we have 



(0), 



(25) 



where C(o) = — /s°^ = /iJ*. Substituting (25) into (24) 

and using then the equalities hs = 0,shs and Is = —ilsls, 
following from (23), we obtain 



^sils + \/phs) = C(o). 

Differentiating (26) we find 

Cls = fllUs - y/phs)/{ls + s/phs), 



(26) 



(27) 



whence it follows that the explosion rate goes through 
the maximum Qj^ in the point where 



(28) 



Combining of (26) and (28) with the "frozen" flux con- 
dition 



/if = hsls = J* 



(29) 



enables one to lock the chain and to readily derive the 
scaling law of concentration explosion. Indeed, from (28) 
and (29) we find 

/,f=p-V4yx, lr=p'/^^/Z. (30) 

Introducing then the scaling function 

C = hs/hf = C/ls 
from (26) and (30) we obtain 



f2, = ^=nfF(C), F(C) = ^ 



where 



M 



C(o)/2C = (M/2)/if. 



(31) 



(32) 
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Integrating (31) with the proviso that C,{T = 0) = 1 we 
easily find 



C - 1/c = 2f7f r 



(33) 



From Eq. (33) it follows that the characteristic time scale 
of explosion is determined by the quantity 7} = l/r?^ 
therefore introducing the reduced time T — T /Tf, we 
finally obtain 



C(T) = T + VT+T2 



whence it follows 



1 



(34) 



(35) 



Two striking features of this result are symmetrical uni- 
versalization |T|~^ ^ of Sl^ beyond the scope of 

interval [—Tf,Tf\ and remarkable symmetry 

Substituting (34) into (19) and using (35) we come to 
the scaling law of growth of the explosion-triggered " an- 
tiflux" Jf-x 

Je. = Je'jiV (36) 

where the amplitude at the point of explosion maximum 
and the scaling function 

j{T) ^aj dec{J - e)S{j - 9)1 Ve. 

Jo 

has the asymptotics 

JiJ) = AajV/^, T » 1. 

The coefficients cm and a j are bound by the relation 
a-Aiaj — I/Vtt therefore by satisfying ^7(0) = 1 we find 
aM = H£!l^ ~ 0.956 and aj = ar^^ ~ 0.590. 

Differentiating (36) we find the singular part of the flux 
relaxation rate in the form 

= -Jex/J. = [Tj']MWiT) (37) 

where the amplitude at the point of explosion maximum 

and the scaling function 

W{T) =cw I de/^[l + (T - eff'^ 
Jo 



has the asymptotics 

WiT) = (37rcv^/8) 171-5/2, _j ^ ^ 

Vr(T) = T»l. 



The coefficients cm and cw are bound by the relation 
cmcw = l/v^ therefore by satisfying 1^(0) = 1 we find 
CM = - 1.046 and cw = ~ 0-539. The 

numerical analysis shows that at T = 0.46205 the scaling 
function W^(T) reaches the maximum 



maxiy(T) = 1.15627... 

whence for the amplitude of catastrophe at the point of 
maximum we find 



maxirji] = (1.15627...)[r7i] 



A/. 



(38) 



Eqs. (30), (32), (34) - (37) give a detailed picture of 
the concentration explosion and annihilation catastrophe 
in the asymptotic limit /C — > oo. Substituting into these 
expressions J^, = AAc and marking the asymptotic val- 
ues of amplitudes with the index (a) in full agreement 
with Q we obtain 



/if(a)=p-i/^VAA:, Zf(a)=pi/4y^^ (gg) 



f^f (a) = (Ai/2)p-i/4^AA; 



(40) 



[Tj%i{a) = (1.43340...)p-5/«A-5/4Ai/4^ (41) 
whence adopting (38) 

max[T7i](a) = (1.65739.. .)p-5/«A-5/4Ai/4_ 

From (17), (18), (37), and (41) it follows that the width 
of flux relaxation rate jump does not depend on p 



\n 



It (X 



A-2/5, 



whereas its amplitude 



maxTji ^ [tj%i oc p-'''^{l ~ pf'^^^'^ 

grows rapidly with diminishing p. So, the less is p, i.e. 
the less L-diffusion restrains the explosion development, 
the more brightly the effect is displayed 

According to (20), one of the necessary conditions for 
the scaling description of passing through the point of 
singularity is the requirement 

VlsP > 1 

which implies that for the medium to be considered as 
semi-infinite in the process of explosion, the explosion 
rate must be much beyond the characteristic rate of the L 
particles diffusion. Combining this requirement with Eqs. 



7 



(34) and (35) and using the equality fJf (a) = y/Kj^-p 
one can easily see that applicability limits of the scaling 
description are described by the inequalities 

< C < Vk:, |T| <C \/^. 

A more rigorous requirement to the quantity /C is im- 
posed by the apparent chain of conditions 



and 



I J. 



Ml 



J* 



li/K,^'^ < 1, 



whence it is to be expected that with a growth of JC 
the amplitudes (39) - (41) ought to be reached mainly 
by the law cx /i//C^/^. One of the remarkable analyt- 
ical advantages of the above given approach is that it 
enables one not only to determine the exact asymptotic 
amplitudes (39)-(41) but, also, to answer the question of 
when and how they are reached. A systematic analysis of 
the crossover to the asymptotics (39)-(41) is given in the 
Appendix. Here I shall focus on the main results. The 
central conclusion of the presented analysis is that fi^f^, 
reaches the asymptotic limit Q,'^\a) much faster than 
fif^^, therefore the point of the explosion maximum is 
defined by precisely by the point of the Hhs maximum. 
According to (A7) 

nffjnfia) = 1 - Bn^i/IC'^^ + On(/C-i/2), 
where i?n ~ 0.0318. With taking into account the equal- 



ity n^s = 



[tj^]m and Eq.(37) it yields 



^Z/^Tia) = 1 + (cm/V2 - Bn)n/IC'/^ + ■■■. 
According to (A8) 

[rJ^U/lTj^Uia) = 1 + Bj^i/K^'^ + Oj{K-^'^) 

where Bj w 0.776. From these expressions we con- 
clude that r^lfg always goes to its asymptotics from below 
whereas [rJ^jAf and always go to their asymptotics 
from above. Remarkably, the coefficient B^ appears to 
be anomalously small so that the contribution of /C~^/^ 
term in the case of Sl^, becomes less than 0.01 already 
at K, > 10^. Below I shall present extensive numerical 
calculations in wide ranges of p, A and /C which demon- 
strate excellent agreement between the numerical data 
and the analytical predictions. 

E. Universalization of annihilation catastrophe. 

Apparently the scaling theory of the annihilation 

catastrophe, developed in the previous section for the 
universal limit no ^ oo, holds in the general case of fi- 
nite no too. Indeed, according to (9) and (10) at finite 
no and A < 1 the chain (24)- (32) remains valid with the 
sole difference that now 



C{o){no) = /iJ*(no)[l - Q{no)] 



become the functions of no. Wc thus have the complete 
scheme to lock the chain (5)-(18) and to answer the ques- 
tion of when and how the universality is reached. It 
remains for us to find the central characteristic of scal- 
ing, namely, the amplitude of explosion fi*^(no), and 
then from Eq.(37) to derive the amplitude of catastro- 
phe {Tj^]M{no)- 

Taking C ^ 1 from Eq.(33) we have 

hs = (/if /2i7f )/|r|. 

Matching this result with Eq.(15) thereby locking the 
chain (5)-(18) we obtain 

= |(1-Q)ftf = |(1-Q)p-V4y:^. 

Using then Eq.(16) we conclude that at A < 1 evolu- 
tion of explosion with growing no is completely defined 
by functions Q(no) and G(no), and find finally the laws 
of universalization of amplitudes of explosion 

flf(no)=f^f(oo)(l + (5a) 

and catastrophe 

[^yVM = [t7VM(i + <5j) 

in the form 

6q = G/2 -Q, 6j = G/4: - 3Q/2. (42) 

Leaving aside details here, I distinguish two main conse- 
quences of (42): 

1) According to (15), (16) the drop of 5a and 5j with 
growing no is surprisingly rapid: 

Q,G ocuQ^ip < Pc), Q,G <x.nQ^{p> pc), 

where 3 < x(p) < 8. Comparing this with the relatively 
slow decrease of ^s[Eq.(13)] 



5s{P < Pc) OC no Saip > Pc) <X Uq 



x/(x-M) 



we conclude that universalization of explosion occurs long 
before /i""" has reached /i^; 

2) According to (42) in the range p < pc with decreas- 
ing p some critical values ^ {p* ) arc reached at which 

6a and 6j reverse their sign ( > +) so that contrary to 

an intuitive reasoning at p < p^^ and p < p*j the ampli- 
tudes fif and maxrj^, respectively, drop with growing 
no- From (15), (16) and (42) we obtain 



and 



?8,n = -Vl7, p*a = 0.0609 



-5/53, p} = 0.0217. 



Note, that this correlates with the behavior of the func- 
tion 6r that pass through zero ( > -|-) at 



J^(no) = AAc[l-hG(no) 



Pr 



1/9. 



8 



IV. NUMERICAL CALCULATIONS 

To test analytical predictions I have carried out ex- 
tensive numerical calculations of Eqs.(l). The numer- 
ical integration of Eqs.(l) was performed by means of 
the implicit discretization scheme of increased accuracy 
with an additional "fictitious" node at the surface [3, 
0. The scheme allowed performing the calculations in 
the system with strong difference in species diffusivities 
with an accuracy down to 10~^%. The space and time 
steps were changed within the ranges 3 x 10~'*-^3 x 10~^ 
and 10^^ 10"^^, respectively, with the number of time 
steps being 10^ ^ 10^. 

In Figs. 2 and 3 are shown the results for A = 10^ 
and p = 0.25, giving detailed picture of formation of the 
universal explosion with growing no. It is seen that in 
accord with Eqs.(42) already at small departures from 
rig(i? — Re) the transient dynamics (7) terminates in ex- 
plosion that remarkably rapidly becomes universal: fur- 
ther growth of no leads to progressing shift of the criti- 
cal point Ti,(no) (Figs. 2a and 3) without changing the 
explosive dynamics in its vicinity but gradually univer- 
salizing the entire self- acceleration trajectory (Fig. 2b). I 
distinguish two important points which characterize the 
universalization process: 

i) In accord with Eqs. (35) and (37) a symmetrical 
"flash" of the explosion rate VLhs {\T\~'^ ^ T'^) (Fig. 
3a) and an accompanying it sharply asymmetrical jmnp of 
the flux relaxation rate rj^ {\T\-^/^ ^ T'^/^) (Fig. 3b) 
form long before the universalization of the corresponding 
amplitudes, shifting self-similarly with growing uq. 

ii) In accord with Eqs. (13) and (17) as no grows, the 
starting point of catastrophe reaches the level coo (Fig. 
3b) long before the starting point of self-acceleration, 
/i™'", reaches the level (Fig. 2). 

In the insets in Fig. 3 are shown the dependences 
^nino) = il|fs(no)/fi|/s(oo) and Ti,{no) plotted on the 
basis of the dependences 57//>;(r, no), represented on the 
main panel and the analogous dependences obtained in 
a wide p range for A ~ 10^ (according to the Appendix, 
here and in what follows the critical point is defined to 
be the point of maximum of fins)- The inset in Fig. 3a 
demonstrates that, in accord with (42), as p decreases, 
the behavior of ri|f^(no) changes qualitatively: a.t p > p^ 
the explosion amplitude r2|/^(no) grows monotonously 
with growing no, reaching f2|f^(cxD) from below whereas 
at p < Pq the explosion amplitude first goes through a 
maximum and then drops with the growing no, reaching 
ilffg{oo) from above. The inset in Fig. 3b demonstrates 
that at all p with the growing uq the numerically calcu- 
lated Ti, values come to the function t"(p, no) calculated 
according to Eq. (14). Moreover, in full agreement with 
Eq. (14) at p > p* the numerical T^, values come to r" 
from below whereas at p < p* the numerical t^, values 
come to r" from above. Below I shall present the results 
of the detailed numerical study of the catastrophe uni- 
versalization in a wide range of p and A and compare 
them with the predictions of (13), (14), and (42) to plot 



on their basis the complete uq — p diagram of universal- 
ization. Before discussing these results, my main goals 
will be: 

1) To demonstrate numerically how with a growth of 
A in the vicinity of the critical point a singularity 
forms, and to show that with a growth of A and no the 
time dependences flnsC^) and [tJ^]{T) collapse to the 
predicted scaling functions S'(T)(35) and M^(T)(37). 

2) Making use of Eqs. (42) for selecting the n^{p) 
region where the contribution of no can be excluded with 
the required accuracy, to study numerically the behavior 
of the r^lfj, and [tJ^]m amplitudes in a wide range of p 
and A, and to demonstrate that with a growth of /C they 
reach their asymptotic values f^s^(a) (40) and [rj^]i^[{a) 
(41) in accord with the predictions of (A7)-(All). 

In Fig. 4 are shown the numerically calculated depen- 
dences r2//s(r) that demonstrate the formation of sin- 
gularity with growing A at no = 4 and p — 0.25. An 
analysis of the given data suggests that as A grows, the 
critical point of the explosion maximum t^;(A) rapidly 
((X A~^) comes to the point of singularity r*(oo), calcu- 
lated according to Eq. (14), so that already at A > 10"* 
the ratio ST^,{A)/T^,{oo) becomes less than 0.001. In the 
inset to Fig. 4 are compared the dependences Qhs{t) and 
^Ls{t) calculated numerically at A = 10*, no = 4 and 
p = 0.25. In accord with Eq. (35), at large fig the curves 
are seen to merge in " synchronous" explosion, asymmet- 
rically coming apart away from the critical point. In Fig. 
5a are shown the data of Fig. 4 for A = 10^ 10^ 10^ 10^ 
replotted in the coordinates fins — T where T = t — t-i,, 
Ti, being the point of the explosion maximum. Here are 
also represented the data of Fig. 3a for no = 1.8, 2.4, 3.1 
(A = 10^, p = 0.25). It is seen that in full agreement with 
Eq. (35) i) the rate of explosion VIhs{T) demonstrates 
the remarkable symmetry —T <-> T and ii) beyond the 
[—Tf,Tf] region, unlimitedly contracting with a growth 
of A, the explosion rate comes to the universal law 1/\T\. 
In the concluding Fig. 5b the data of Fig. 5a are rep- 
resented in the scaling coordinates ^Ihs/^hs ~ ^ where 
T = T /Tf = TVl.^^. It is seen that the numerically cal- 
culated dependences are perfectly collapse to the scaling 
function ^(T) (35). 

Let us now turn to an analysis of the flux relaxation 
rate. Fig. 6a demonstrates the dependences rJ'^(T) cal- 
culated numerically for the same parameters as in Fig. 
5a (the data are shown only for no = 4). The data anal- 
ysis suggests that in accord with Eqs. (17), (18), and 
(37) with growing A in the vicinity of the critical point r^, 
forms a singular jump of rj^ the width of which contracts 
unlimitedly by the law \Tcat \ oc A~^/^ and the amplitude 
of which grows unlimitedly by the law maxTj ^ cx A^/'^ 
(see below). Note that in accord with Eq. (37), after 
the critical point has been passed (T ^ 7}), the relax- 
ation rate drops by the A independent law cx 1/ \JTp^ 
reaching at times T ~ p/wq L-diffusion-controUed 
limit wo/p so that in the limit of small p there arises a 
most dramatic consequence of the annihilation catastro- 
phe: an abrupt, practically instantaneous (on the scale of 
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ujq) disappearance of the flux [3]. Based on the data of 
Fig. 6a in accord with Eq. (17) the time dependences 
of the singular part of the flux relaxation rate were cal- 
culated [Ty^](T) = tJ^(T) — a;o(l + w) which were then 
replotted in the scaling coordinates [•''7^](T)/[tJ"^](0)— T. 
The results are demonstrated in Fig. 6b. It is seen that 
in perfect agreement with Eq. (37) with growing A the 
numerical results collapse to the scaling function M^(T). 
For a more detailed illustration in Figs. 6c and 6d the 
data of Fig. 6b are represented in double logarithmic 
coordinates in a wider range of |T| separately for T < 
(Fig. 6c) and T > (Fig. 6d). Ibid are given the data 
based on the numerical calculation for Uq = 1.8,2.4,3.1 
(A = 10^, p = 0.25) which demonstrate a collapse to 
the scaling function M^(T) with growing uq. It is seen 
from Fig. 6c that at A = 10^ the range of the scaling 
growth regime of [tJ^] reaches six orders of magnitude. 
As in this case at the starting point of growth the ra- 
tio [tJ^]/ujq ~ 10""*, it implies that the accuracy of the 
numerical calculation of tJ^ reaches 10~^%. Extensive 
numerical calculations in a wide range 10~^ < p < 1, a 
part of which will be given below, have shown that at all 
the investigated p values with a growth of A (and there- 
fore JC) the obtained numerically normalized dependences 
^HsC^) and [t7^](T) collapse, respectively, to the scal- 
ing functions S{T) and VF(T). We thus conclude that 
the scaling theory of catastrophe perfectly agrees with 
the numerical results. 

Let us now come to a numerical study of the behav- 
ior of the amplitudes fl^g{p, A, no) and [tJ^]m{P: A, no) 
which are the central characteristics of the scaling regime 
of catastrophe. Following the above stated program, I 
shall begin with the results derived in the universal limit 
uq ~^ oo. The numerical calculations were performed 
within A = lO-^^lO^ for p = 0.01, 0.03, 0.1, 0.25, 0.5, 0.75. 
In order that the contribution of the initial conditions 
be excluded, the initial number of particles no was de- 
pending on p selected from the range no = 10 200 
so that in accord with Eqs. (42) this contribution may 
not exceed 10~^%. In Fig. 7a are shown the numer- 
ically calculated dependences 7^ = Sl^g/f2*^(a) and 
7.7 = [tJ^]m /[Tj^]M{a) as functions of /C//i^. It is seen 
that with growing K. the numerically calculated ampli- 
tudes and [tJ^Jm come, respectively, to the asymp- 
totic values ri^(a) and [Tj^]M{a) calculated analytically 
according to Eqs. (40), (41). Remarkably, in accord with 
the predictions of Eqs. (A7), (A8) i) 7n comes to 1 from 
below whereas ■jj comes to 1 from above; ii) 70 comes to 1 
much faster than 7,7; iii) the law by which 77 approaches 
1 in a wide range of K, //i* is described with excellent ac- 
curacy by the principal term of Eq. (AS). For a more de- 
tailed illustration of iii) in Fig. 7b in double logarithmic 
coordinates are presented the dependences ej = 77 — 1 vs. 
K,/ jji^. Ibid are shown the numerically calculated depen- 
dences = {Ji, — Jm)I J* vs. /C//i^ which demonstrate 
the law by which Jm approaches J* with growing /C. It is 
seen that the numerically calculated ej and values at 
p not too close to 1 perfectly fall on the analytic depen- 



dences BjixlK}!^ [Eq. (AS)] and B^iilK}!^ [(Eq. (A4)], 
respectively, shown in thick lines (note that according to 
(AlO) at p — > 1, /i ~ 1 —p the dominant term in ej to 
/C oc /i"* ^ 00 becomes cx l/v^)- In accord with Eqs. 
(A7), (A9), (All) owing to the anomalous smallness of 
the coefficient the value of = 1 — 70 ought to de- 
crease with growing /C by the law cx p/^/lC or (at small p) 
even faster down to very low values of eo ~ 10~^. To il- 
lustrate these predictions in Fig. 7c are shown the depen- 
dences of en on ^/K. in double logarithmic coordinates. 
At comparatively high p the numerically calculated eo 
values are seen to fairly fall on the analytic dependences 
shown in dashed lines. As p decreases the effective "rate" 
of drop grows only insigniflcantly, the en value itself 
becomes very small already at K. ^ 10^ — 10^, therefore 
the analytical description of en in this region necessitates 
additional terms. Summarizing, we conclude that in the 
universal limit the represented theory gives an exhaustive 
picture of evolution of catastrophe and explosion ampli- 
tudes. 

It only remains for me to complete this section by 
demonstrating the results of extensive numerical study 
of the regularities of catastrophe universalization with 
growing no. I have studied the behavior of the de- 
pendences T^(no), r2|fg(no ),[t7VK) and /i™(no) at 
"scanning" no from n§ to 10^ in wide ranges of A = 
10^ 10*^ and p = lO'^ 0.97. Based on the obtained 
data for each of the studied p and A values I calcu- 
lated the dependences drijiQ) = Ti,{nQ)/T^ — l,(5si("-o) = 
ra(no) - l,Sj{no) = [rj^lif ("-o)/[t7^]m(oo) - 1 and 
Ss{no) — ft.™™ (no) /Ac — 1 which I then compared with 
the analytic predictions. I have found that in the region 
of small (5i(no)(i — r, £7, J, s) the behavior of functions 
(5i(no) is described with remarkable exactness by Eqs. 
(14), (42) and (13). Fig. S represents the concluding 
no — p diagram of universalization where are compared 
the positions of the boundaries \5i\ = 0.01(i = r, fi, J) 
and Ss = 0.1 resultant from a great number of numerical 
data (some of which are given in the inset) for A — 10^ 
(in the case i = J for A = lO"^) and calculated from (42), 
(14) and (12). Excellent agreement of the analytic and 
numerical results {not shifting with the further growing 
A [9]) needs no comments. 



V. DISCUSSION AND CONCLUSION 

Finite-time singularities - blowup solutions develop- 
ing from a smooth initial conditions at a particular time 
- provide probably the most dramatic manifestation of 
strongly nonlinear effects that can occur in nature |lOj |. 
The formation of finite-time singularities is observable in 
a wide spectrum of nonlinear systems (Jang-Mills fields 
11 1 , black holes 0| , self-gravitating Brownian particles 
13 1, turbulent flows fl3. iet eruption chemotaxies 
16 1 , and earthquakes [13 to name only a few) therefore 
the description of scenarios of finite-time singularities de- 
velopment is a fundamental problem which attracts a 
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wide interdisciplinary interest. 

In this paper for the first time a systematic theory of 
formation of the universal annihilation catastrophe from 
a smooth initial distribution has been developed and ex- 
tensive numerical calculations of the regularities of the 
catastrophe formation in a wide range of parameters have 
been presented. The main results may be formulated as 
follows: 

1) The exact condition of the i/-diffusion-controlled 
annihilation regime has been found; 

2) The rigorous proof of the threshold arising of self- 
acceleration has been given; 

3) The exact expression for the time point of catastro- 
phe has been obtained; 

4) The closed scaling theory of passing through the 
point of singularity has been given and the scaling laws 
of concentration explosion and annihilation catastrophe 
have been derived; 

5) The laws of universalization of concentration explo- 
sion and annihilation catastrophe have been found; their 
surprisingly rich "structure" has been revealed; 

6) A remarkable agreement has been found between 
the analytical predictions and the results of numerical 
calculations. 

Summarizing, I believe that the analysis presented may 
pretend to be one of the most striking examples of de- 
tailed description of formation of finite-time singularity. 
I shall distinguish here two most bright features of anni- 
hilation catastrophe: 

i) In the majority of the models which demonstrate 
the formation of finite-time singularities an analytical de- 
scription of the singularity development (based on prop- 
erties of self- similarity) appears possible only to some 
narrow vicinity of the critical point beyond which the 
solution cannot as a rule be continued or is principally 
impossible. One of the main advantages of the here 
presented theory is the asymptotically exact scaling de- 
scription of passing through the point of singularity which 
yields a complete dynamical picture at the both sides of 
the critical point; 

ii) Arising as a result of explosive growth of the "an- 
tiflux" Jex at the background of slow relaxation of the 
diffusion-controlled flux the annihilation catastro- 
phe demonstrates a peculiar singular behavior at which 
two explosive processes {Qhs and flLs) are developing 
simultaneously, effectively "compensating" one another 
so that for an external observer of fiux (J) the explosion 
dynamics goes unnoticed up to the critical point r^, in 
the vicinity of which "decompensation" of explosions is 
manifested as a sudden singular jump of the flux relax- 
ation rate. In the limit of small p this brings about a 
most radical consequence - an abrupt disappearance of 
the fiux. 

Let us discuss in conclusion the conditions and pos- 
sibilities of an experimental observation of the annihila- 
tion catastrophe. The irreversible bimolecular reaction 
^ -|- B ^ is one of the most abundant reactions there- 
fore it is to be expected that the predicted phenomena 



can, in principle, be observed in a wide class of physi- 
cal, chemical and biological systems with a "catalytic" 
interface which, because of the high energetic barrier, 
does not let diffusing particles A go from medium 1 to 
medium 2 and diffusing particles B from medium 2 to 
medium 1 so that the reaction A + B can occur only 
at the interface between the media |18|| . p^ . Leaving 
aside here discussing systems of such type, I shall focus 
on the main object of the model in question, namely, ad- 
sorption - desorption systems (Fig. 1). Until now most 
of the theoretical studies on the Aads + Bads cat- 
alytic reaction (the Langmuir-Hinshelwood process which 
is also often referred to as the monomer-monomer cat- 
alytic scheme) have been performed under the assump- 
tion that diffusion into the bulk can be neglected (w| 
- 28] and references therein). Such an assumption is 
valid in low-temperature systems with high surface-bulk 
crossover barriers, i. e., in systems with the negligibly 
small bulk solubility of A and B particles. Here I address 
to the wide class of catalytic systems where the surface- 
bulk crossover barriers are not too high and, therefore, 
adsorption-desorption processes are always followed by 
a more or less intensive diffusion of A and B particles 
into or from the bulk where reaction between A's and 
B's is energetically forbidden :29]. This class of catalytic 
systems is not only of fundamental interest for surface 
science but, also, of a considerable applied interest for 
describing the interaction kinetics of gases with metals 
at high temperatures {^2^ - [s^l and references therein). 
In the work 33] the theory has been developed for the 
diffusion-controlled associative desorption of like parti- 
cles, Aads + Aads -> A2.gas ^ 0, from dissolved state 
into vacuum. Adopting this theory, a complete picture 
of diffusion-controlled thermodcsorption of hydrogen and 
nitrogen has been constructed in good agreement with 
the available experimental data. The theory reported in 
this work gives a systematic description of the diffusion- 
controlled kinetics of associative desorption into vacuum 
of unlike particles 

Aads + Bads ^ ABgas ^ 0, 

which are initially uniformly dissolved in the bulk. I shall 
focus here on discussing a possibility of observation of the 
predicted effects for one of the most important surface re- 
actions of carbon monooxide CO thermodcsorption from 
metals into vacuum 

Cads + Oads *■ COgas * 0. 

It is to be mentioned first that the continual description 
(1) holds as long as the "diffusion length" of the explo- 
sion at the point of maximum remains much greater than 
the monolayer thickness a Q, Sxm '3> a/£, 

whence there follow the limitations 

r!f«(^/a)2, /C«p2(^/a)4. 

Taking, for example, £/a ^ 10^ and p ~ 0.01 we come to 
the requirements il*^ ^ 10^ and K, ^ 10® to see that at 
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any value of the reaction rate constant k the specimens 
must have macroscopic sizes in order a considerable ef- 
fect be observed. Based on the data of monograph [3?il |. 
I shall make estimations for three refractory metals, i.e. 
niobium, tantalum, and molybdenum which at elevated 
temperatures dissolve carbon and oxygen in quite large 
amounts. According to 30] , at temperatures of intensive 
thermodesorption of CO in the range from T ~ 1600°C 
to melting point for coefficients of carbon and oxygen 
diffusion in these metals we find, respectively, Dq ~ 



)crn? /s, 
-1 



(10-^ 10-5)ctoVs and Dq (10"^ 10" 
whence it follows p = DcJDq ^ 10"^ ^ 10"^ Ac- 
cording to the data of [33 - |32| the desorption rate 
constant of CO in said temperature range alters within 
K ^ {lQ~^^^10~^^)cm'^ / s. Substituting these values into 
the expression 

A = 6c{0)ni/Dc 

and taking Sc{0) = cc(0) - co(0) ~ W°cm-^ and 
H. ~ 0.1cm we find that in said temperature range the 
A parameter value changes within A ~ 10^ 10^. For 
the density of the diffusion-controlled desorption flux of 
CO at the critical point we find h ^ DcSc{0)/£ ^ 
10^^ 10^^ particles /cm^s. We thus conclude that in 
a study of isothermal desorption of CO at elevated tem- 
peratures under high vacuum the predicted sharp jump 
of the flux relaxation rate can confidently be registered 
experimentally with a standard measuring technique. 
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APPENDIX A: CROSSOVER TO THE SCALING 
REGIME OF CATASTROPHE. 

From (40) and (41) it follows that in the limit of large 
IC oo the ratio 



1/4 



0. 



nf 

I shall show below that reaches the asymptotic limit 
(40) much more rapidly than therefore it is the 

point of maximum of Qhs that defines the point of max- 
imum of the explosion. With allowance for the fact that 
^Ls/^Hs 1 only asymptotically, at finite K. instead of 
(26) one has to write 



-{0)- 



(Al) 



whence at the point of explosion maximum fins = we 
find 



where for completeness the term with the derivative is 
held 



C. 



(0) 



and it is taken into account that from the condition 
^Hs = it follows (Ils = [tJ^]. Differentiating (37) 
and calculating the arising integral, we find the leading 
term in [tj in the form 

where jm = 3r(3/4)V27r ^ (3/4)aM ~ 0.717. Sub- 
stituting this result into Eq. (A2), using the equality 
flf^ = njf^ +ljo + [tJ^]m, and taking « Of, after 
separating out the leading terms we find 



^hf/lf = 1 + Brii/IC^'^ + Or{IC-^'^) (A3) 



where 

^ 2cM ~1M ^ Q g^2^ 
V2 

At the point of the explosion maximum allowing for the 
contribution of J^x (36) instead of (29) one has to write 

Substituting (36) here we have 

hflf - J.[l - B^^i/K.^/^ + 0.(/C-i/2)]^ (A4) 

where 

= aM/\/2« 0.676. 
Eqs.(A3) and (A4) immediately give 

hf/hf{a) = 1 + Bhti/JC^'^ + OhiJC-^'^). (A5) 

where 

/if(a)=p-i/4v^ 

and 

{2cM - aM - 1m) 



Bh 



2V2 



0.148, 



and 



Is' /If (a) = 1 - i3iM//Ci/4 + Oz(/C-i/2), (A6) 



where 



and 



;f(a) = pi/4yx 
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From (A5) and (A6) it follows that ft,*^ always comes 
to its asymptotics (a) from above whereas if^ always 
comes to its asymptotics l^^ (a) from below. Essentially, 
the coefficient Bh is much smaller than Bi and, therefore, 
with growing /C the asymptotics h^{a) is reached much 
earlier than the asymptotics l^{a). Substituting then 
(A3) into (Al) and using (A5), (A6) we obtain 



n^jnfia) = 1 - Bnii/K}'^ + On(/C-V2), (A7) 



where 



and 



2s/2 

From (A7) it follows that fl^^ always comes to its asymp- 
totics 0^(a) from below. Remarkably, the coefficient Bq 
appears so small that already at /C > 10^ the contribution 
of the K~^/^ term becomes less than 0.01. Substituting 
now (A4), (A5), and (A7) into the expression 

for the amplitude of the catastrophe at the point of ex- 
plosion maximum we find 

[t7 V/[^7 = 1 + Bjti/K}'^ + Oj{K-^l^),{k%) 

where 

[^7V(«) = (0.369834...)/.3/2p-5/8ji/4 



and 



Bj = 



au + (1/8) (2cm - Qm) 
^/2 



0.776 



From (AS) it follows that [tJ^]m always comes to its 

asymptotics [tj ]m{o) from above and, due to the com- 
paratively high Bj value, reaches its asymptotics much 
slower than . 

The expressions (A5), (A6), (A7), and (AS) completely 
define the principal picture of the crossover to the scal- 
ing catastrophe and explosion regime in the limit of large 
/C — ^ oo. To be complete, I shall calculate now the correc- 
tions Oa{K~^/'^) and Oj{IC^^^^) which, as one can easily 
see, may be of two types 0\{p/^/lC) and Of (/U^/\/^)- As 
/i~l— p^Oatp— *1 the corrections appear essen- 
tial at p close to 1. The calculations give 



OP = -(7rV4)p/V^ w -2.467p/\/^ 



(A9) 



and 



(4aM/cM - ^ +0.S09p/v^. (AlO) 



Due to the anomalously low value of the coefficient 
I shall also give the term 0[^. The calculations yield 
= [Bn(3Bn/2-Sft-B*/S)-B2/128]MVV^ whence 
after substituting the coefficients we find 



-0.00946/i7v^- 



(All) 
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